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This  document  presents  a discussion  of  the  various  types  of  wind  shear 
which  cause  problems  to  aircraft  operations.  The  presentation  is 
limited  to  lew-level  wind  shear  such  as  could  be  encountered  in  the 
vicinity  of  airports  by  aircraft  on  approach  to  landing  and  on  takeoff. 
The  types  of  shear  discussed  are  primarily  limited  to  frontal, 
thunderstorm  and  shears  associated  with  stable  and  neutral  boundary 
layers.  The  discussions  emphasize  modeling  of  the  wind  shear  for 
simulation  purposes. 

Wind  shear  is  discussed  not  only  as  a change  in  absolute  wind  speed 
but  ail  so  as  a result  of  wind  direction  change.  Same  insight  is 
given  relative  to  a combination  of  wind  speed  and  direction  change. 

It  is  to  be  noted  that  the  condition  affecting  aircraft  operations  is 
not  one  shear  parameter  but  is  a combination  of  several,  e.g. 
horizontal  shear,  vertical  shear,  wind  direction  change  and 
height  of  shear  above  ground  level. ^ 
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1.0  INFLUENCE  OF  WIND  SHEAR  ON  AIRCRAFT  OPERATIONS 
1.1  Introduction 

Wind  is  an  important  consideration  in  the  analysis 
of  airplane  flight  in  the  atmospheric  boundary  layer,  both 
because  of  short  scale  gusts  or  turbulence  and  because  of  large 
scale  variations  of  the  mean  wind.  In  the  planetary 
boundary  layer,  the  mean  wind  decays  toward  the  ground 
and  has  considerable  horizontal  variations  due  to 
irregularities  in  terrain.  Thus,  both  spatial  and 
temporal  variations  occur  in  mean  surface  winds  encountered 
along  ascending  and  descending  flight  paths. 

Many  analyses  of  airplane  motion  have  been  carried 
out  (1-1,  2,  3)  which,  in  general,  consider  only  constant 
winds  and  thus  neglect  effect  of  wind  shear.  Recent 
reports  of  commercial  airline  accidents,  however,  have 
identified  wind  shear  in  the  planetary  boundary  layer 
as  being  a significant  factor  in  the  cause  of  these 
accidents,  illustrating  that  variable  winds  are  an 
imports  design  consideration. 

This  report  begins  with  a review  of  the  reported 
influence  of  wind  shear  on  aircraft  accidents.  Recent 
major  accidents  where  wind  shear  was  a primary  factor  are 
described . 
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State  of  the  art  mathematical  models,  tabulated  data  of  mean 
wind  profiles  and  fields  which  can  be  generally  used  in 
computer  simulation  of  aircraft  flight  through  the  wind 
shears  associated  with  these  meteorological  conditions  are 
presented. 


1 . 2 Problems  of  Wind  Shear  in  Aircraft  Operations 

Accidents  caused  by  encountering  strong  wind  shear, 


turbulence  and  gusts  during  terminal  flight  operations  are 
becoming  increasingly  more  evident  as  improved  documentation 
o.f  these  conditions  is  achieved  through  advanced  instrument 
technology,  such  as  digital  flight  data  recorders.  Reference 
[1-4]  reports  that  weather  conditions  were  significant 
factors  in  air  carrier  accidents  within  the  last  two  years 
and  Wyatt  [1-5]  reports  weather  to  be  a factor  in  approxi- 
mately 39%  of  all  fatal  accidents  which  occurred  in  general 
aviation  between  1964  to  1973. 

To  what  degree  these  were  due  to  wind  effects  is  un- 
certain. However,  many  missed  approaches  which  have  simply 
been  classified  as  pilot  error  may,  in  the  light  of  today's 
knowledge,  be  traceable  to  wind  shear. 

Recently,  a number  of  well-defined  incidents  where  wind 
shear  has  been  a major  factor  contributing  to  serious 
accidents  have  been  documented.  Laynor  [1-6]  reports  that 
an  Iberian  Airlines'  DC  10-30  crashed  on  December  17,  1973, 
about  500  feet  short  of  the  runway  while  making  an  Instrument 
Landing  System  (ILS)  approach  to  Logan  International 
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Airport,  Boston,  Massachusetts. 

The  meteorological  data  indicated  that  the  winds  aloft 
at  altitudes  above  1000  feet  in  the  Boston  area  were  generally 
from  the  south  at  speeds  approximately  40  kts.  The  surface 
wind  measured  at  Logan  Airport  was  from  a westerly  direction 
at  9 kts.  These  data  suggest  that  wind  shear  existed  at  the 
time  of  the  accident. 

Examination  of  the  data  recorded  on  the  aircraft  digital 
flight  data  recorder  generally  verified  that  the  aircraft 
penetrated  an  altitude  band  where  a sudden  change  in  wind 
direction  and  speed  occurred.  The  results  of  re cons true ting 
the  approach  in  a Douglas  DC-10  simulator  disclosed  that  a 
wind  shear,  characterized  by  a diminishing  tailwind  component, 
encountered  at  low  altitudes  during  the  period  of  approach 
where  the  pilot  is  transitioning  from  automatic  to  manual 
flight,  produced  a situation  in  which  the  aircraft  descended 
below  the  glide  path.  Winds  representative  of  those 
experienced  during  the  accident  are  shown  in  Table  1-1. 

On  January  4,  1971,  a FAA  Douglas  DC-3C,  N7  crashed  at 
LaGuardia  Airport,  New  York,  approximately  2000  feet  short 
of  the  approach  threshold.  The  probable  cause  of  the 
accident  was  reported  [1-7]  as  the  failure  of  the  pilot  to 
recognize  the  wind  shear  conditions  and  compensate  for  them. 
The  meteorological  conditions  which  prevailed  the  day  of  the 
crash  were  a warm  front  between  the  Kennedy  and  LaGuardia 
airports  and  the  winds  aloft  were  southwesterly  and  quite 
strong  while  the  surface  winds  were  from  the  northeast  and 
relatively  gentle. 
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TABLE  1-1 

REPRESENTATIVE  WINDS  ALONG  FLIGHT  PATH  OF  DC  10-30  AIRLINER 
WHICH  CRASHED  AT  LOGAN  AIRPORT,  DECEMBER  19,  1973  [1-6] 


Altitude 

ft. 

Longitudinal 

Component 

kt 

Lateral 

Component 

kt 

Direction/Speed 

North/kt 

1600 

25  tail 

26  left 

169/36 

1000 

23  tail 

26  left 

176/35 

650 

21.5  tail 

25  left 

176/33 

500 

18.0  tail 

23  left 

179/29 

420 

12.5  tail 

18.5  left 

183/23 

350 

9 tail 

14.5  left 

190/17 

255 

3 tail 

10  left 

229/9 

180 

6 head 

2 left 

270/5 

100 

6 head 

2 left 

296/6 

0 

4 head 

2 left 

296/5 

Sowa 

[1-8]  reports  four 

cases  of  wind  shear  ranging 

from  a catastrophic  accident 

to  a frightening 

experience 

for  which 

the  primary  cause 

can  be  attributed 

to  wind  shear. 

Figure  1-1  from  [1-8]  illustrates  the  flight  recorder  traces 
of  two  dissimilar  aircrafts  at  widely  separated  regions  of 
the  earth.  The  figures  suggest  that  the  aircraft  most  likely 
encountered  a strong  headwind  near  3000  feet  resulting  in  a 
pitched-up  attitude  and  a high  rate  of  ascent.  At  8000  feet 
with  the  throttle  and  altitude  probably  adju. „ed  to  trim  out 
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the  high  climb  rate,  a tailwind  was  presumably  encountered, 
resulting  in  a very  low  pitch  attitude  and  a rapid  rate  of 
descent  over  a distance  on  the  order  of  7000  feet.  Another 
1000  feet  and  the  incident  may  well  have  been  fatal. 

Nine  missed  approaches  (see  Table  1-2  from  [1-8] ) made 
at  J.  F.  Kennedy  Airport,  New  York,  between  2152Z  and  2354Z 
on  January  4,  1971,  highlight  a problem  also  associated  with 
wind  shear.  Chambers  [1-9]  reports  a warm  front  approached 
JFK  from  the  south  arriving  approximately  2300Z.  At  2330Z, 
winds  were  as  indicated  in  Table  1-3. 

The  wind  speeds  shown  demonstrate  that  a wind  shear  of 
2 kts/100  ft  was  present.  Prior  to  the  passage  of  the  front, 
the  surface  wind  was  040°/7  kts.  During  approach  under  these 
shear  conditions,  a decreasing  tailwind  with  decreasing 
altitude  was  encountered  and  heavy  (slow  response)  aircraft 
were  unable  to  prevent  a speed  overrun.  They  were  also  high 
on  the  glide  slope.  This  conceivably  caused  many  of  the 
missed  approaches  shown  in  Table  1-2.  Difficulties  occurred 
during  takeoff  as  well  since  departing  airplanes  experienced 
loss  of  air  speed  due  to  the  wind  shear  and  also  a thrust 
loss  due  to  flying  into  warmer  air  aloft. 

Chambers  [1-9]  points  out  that  British  Overseas  Airline 
Company  pilots  have  reported  other  cases  of  marked  wind  shear 
at  New  York.  Additionally,  he  notes  that  BOAC  lost  an  airplane 
shortly  after  takeoff  from  Kano,  Nigeria  in  1965  due  to  wind 
shear  from  a cold  air  down  draft  accompanying  a thunderstorm. 
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FIGURE  1-1  FLIGHT  RECORDER  TRACES  OF  TWO  DIFFERENT  AIRCRAFT 
PASSING  THROUGH  UNEXPECTED  WIND  SHEAR  U-8], 


INDICATED  AIRSPEED  (KTS)  INDICATED  AIRSPEED  (KTS) 


• TABLE  1-2 

DATA  FROM  TOWER  LOG  OF  A MAJOR  AIRPORT  [1-8] 


Time  ( Z ) 

Aircraft 

Comment 

2152 

Twin  Turbo 

Missed 

approach  & diverted 

2200 

Wide-body 

Landed 

2237 

4-Engine  Jet 

Missed 

approach 

‘ 2300 

4-Engine  Jet 

Landed 

second  approach 

2302 

Wide-body 

Missed 

approach  & diverted 

2304 

4-Engine  Jet 

Missed 

approach 

2329 

Tri-Jet 

Missed 

approach 

2333 

4-Engine  Jet 

Missed 

approach 

2341 

4-Engine  Jet 

Landed 

second  approach 

2346 

Tri-Jet 

Landed 

second  approach 

2349 

Wide-body 

Missed 

approach 

2353 

4-Engine  Jet 

Landed 

second  approach 

2354 

Wide-body 

Missed  approach 

0013 

Changed 

. from  runway  04R 

to  22 

0020 

Wide-body 

Landed 

second  approach 

0023 

Wide-body 

Landed 

second  approach 

TABLE  1-3 

WINDS  OVER  JFK  ON 

JANUARY  4,  1971  [1-9] 

Direction/Speed 

(kts) 

Altitude 

220/05 

Surface 

215/26 

1000  ft. 

220/46 

2600  ft. 

7 

1 . 3 Effects  of  Wind  Shear  on  Aircraft  Flight 

Basically,  there  are  two  potentially  hazardous  shear 
situations  [1-10] . First,  a tailwind  may  shear  to  either 
a calm  or  headwind  component.  In  this  instance,  initially 
the  airspeed  increases,  the  aircraft  pitches  up  and  the 
altitude  increases  (Figure  1-2).  Second,  a headwind  may 
shear  to  a calm  or  tailwind  component.  Initially,  the 
airspeed  decreases,  the  aircraft  pitches  down  and  the  altitude 
decreases  (Figure  1-3).  Aircraft  speed,  aerodynamic  char- 
acteristics, power/weight  ratio,  powerplant  response  time  and 
pilot  reactions  along  with  other  factors  have  a bearing  on 
the  severity  of  wind  shear  effects. 

The  potential  hazards  of  wind  shear  therefore  suggest 
that  manned  flight  simulators  should  be  programmed  to  train 
flight  crews  to  cope  with  shear  conditions  and  for  fast  time 
computer  analysis  to  relate  the  potential  hazards  posed  by 
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the  wind  shear  phenomenon  to  different  types  of  aircraft  and 
control  systems.  Section  2.0,  which  follows,  identifies  the 
various  wind  shear  regimes  while  Sections  3.0,  4.0  and  5.0 
provide  a concise  review  of  the  state-of-the-art  knowledge  of 
frontal  wind  shears,  thunderstorm  wind  shears  and  wind  shears 
found  in  stable  and  neutral  atmospheric  boundary  layers,  res- 
pectively. Formulation  of  existing  models  and/or  wind  data 
into  mathematical  expressions  is  given  in  these  sections  for 
later  applications  to  the  engineering  models  to  be  developed 
in  the  final  report. 
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2.0  TYPES  OF  WIND  SHEAR 
2 . 1 Introduction 

Wind  shear  is  caused  by  several  different  motions'  of 
atmospheric  air  masses.  The  atmospheric  boundary  layer  in 
its  natural  state  always  contains  some  degree  of  wind 
shear.  The  geostrophic  wind*  aloft  which  flows  parallel  to 
the  isobars  is  turned  and  retarded  by  frictional  forces  as 
the  earth  surface  is  approached.  This  creates  the  atmospheric 
turning  layer  which  has  both  directional  and  speed  variation 
with  height  resulting  in  vertical  wind  shear.  Near  the  earth, 
horizontal  variation  in  winds  occur  due  to  terrain  irregu- 
larities, creating  horizontal  wind  shear.  Thus,  during  most 
approach  and  takeoff  operations,  some  degree  of  wind  shear 
is  encountered.  The  strength  of  the  shear  and  the  degree  to 
which  it  becomes  hazardous  is  dependent  upon  the  existing 
combination  of  meteorological  conditions.  Under  conditions 
of  strong  nighttime  temperature  inversions,  which  may  have 
in  addition  a low  level  jet  near  the  top  of  the  inversion, 
severe  shears  can  be  created,  whereas,  under  conditions 
of  a super  adiabatic  lapse  rate,  shear  may  be  destroyed  by 
vigorous  turbulent  mixing.  Wind  shear  can  occur  in  land/sea 

*The  qeostrophic  wind  results  from  a balance  of  the  horizon 
tal  pressure  gradient  and  the  Coriolis  force,  and  it  blows  along 
straight  parallel  isobars  above  the  boundary  layer. 
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breezes  and  in  anabatic/katabatic  winds  associated  with 
local  topography. 

Large  wind  shear  occurs  in  association  with  cold  air 
downdrafts  spreading  over  the  ground  outwards  from  a thunder- 
storm (gust  front)  which  may  precede  the  storm  by  as  much  as 
10  miles.  Slowly  moving  warm  and  cold  fronts  create  wind 
shear  primarily  of  a directional  variation  but  also  at  times 
of  a wind  speed  variation.  Cases  of  low  level  jets  with 
accompanying  strong  shear  have  been  noted  in  warm  sectors 
near  cold  fronts. 

Sections  3.0,  4.0  and  5.0  of  this  report  provide  in- 
sight into  the  current  understanding  and  prediction  schemes 
available  for  the  three  major  sources  of  wind  shear,  i.e., 
frontal  wind  shear,  thunderstorm  wind  shear  and  wind  shear 
associated  with  stable  and  neutral  boundary  layers. 


- S*?i  JgSfifeStoa*- . 


3.0  FRONTAL  WIND  SHEAR 


Fronts, as  described  in  this  section,  refer  to  warm  and 
cold  frontal  surfaces  in  motion  and  not  to  the  thunderstorm 
gust  front  discussed  later.  Not  all  fronts  produce  signifi- 
cant shear  and,  in  fact,  many  have  gradual  transitions  in 
wind  speed  over  broad  transition  regions.  Those  cold  and 
warm  fronts  which  do  have  sharp,  narrow  transition  zones 
generate  significant  wind  shear  that  is  particularly  dis- 
turbing due  to  associated  bad  weather.  It  was  wind  shear 
created  by  a slowly  moving  warm  front  that  caused  the 
numerous  missed  approaches  reported  at  J.  F.  Kennedy  Airpoit 
January  4,  1971. 

Mathematical  expressions  for  the  structure  of  wind  fields 
accompanying  fronts  are  not  presently  available,  but  criteria 
for  those  fronts  which  contain  significant  wind  shear  is  given 
by  Sowa  [1-8]  as:  1)  fronts  which  have  a 10°F  (5°C)  or  more 

temperature  difference  immediately  across  the  front  at  sur- 
face level;  and/or  2)  fronts  which  are  moving  at  30  kts  or 
faster. 

Wind  shear  associated  with  fronts  can  be  either  direc- 
tional shear  or  speed  shear.  Figure  3-1  illustrates  two 
cases  of  directional  shear.  The  large  solid  line  in  Figure 
3-la  indicates  the  cold  front.  The  wind  below  the  front  at 
A (surface  wind)  is  from  320°  and  the  wind  above  the  front 
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(which  is  paralled  to  the  isobars)  is  from  220®.  Similarly, 
point  D in  Figure  3-lb  which  has  a warm  front  sloping  over 
it,  has  a surface  wind  from  040°.  Above  the  front,  the  wind 
is  paralled  to  the  isobars  from  230°.  The  height  above  the 
landing  surface  at  which  the  directional  change  occurs  is 
difficult  to  determine.  A cold  front  moving  at  30  kts  or 
more  might  reasonably  be  expected  to  have  a frontal  slope 
such  that  the  front  is  5000  feet  above  the  airport  three 
hours  after  the  passage  of  the  front.  Figure  3-2  shows  a 
linear  extrapolation  of  the  above  observation.  The  shear 
associated  with  a warm  front  usually  exists  below  5000  feet 
and  for  approximately  six  hours  prior  to  the  warm  front 
passing  the  airport.  The  shear  becomes  negligible  once  the 
warm  front  has  gone  by. 

Wind  speed  changes  as  opposed  to  wind  directional 
changes  occur  most  frequently  with  slow  moving  warm  fronts 
having  large  temperature  differences.  Very  little  infor- 
mation on  predicting  these  speed  changes  is  available. 

The  presence  of  turbulence  with  the  frontal  wind  shdar 
may  or  may  not  occur.  Figure  3-3  indicates  the  approximate 
conditions  of  speed  and  temperature  difference  which  will 
induce  turbulence  along  with  the  wind  shear. 

Recent  data  available  from  the  Atmospheric  Variability 
Experiment  Number  IV  (AVE  IV)  provides  three  hour  soundings 
of  wind  speed  and  direction  at  25  m increments  in  elevation. 
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TEMPERATURE  DIFFERENCE  WHICH  WILL  INDUCE 
TURBULENCE  ALONG  WITH  WIND  SHEAR/  I1_83* 
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These  data  are  available  for  23  stations  throughout  the 
eastern  half  of  the  United  States.  Analysis  of  these  data 
to  extract  detailed  structure  of  fronts  and  to  develop  a 
predictive  model  of  frontal  winds  should  be  carried  out. 
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4.0  THUNDERSTORM  GUST  FRONTS 


4.1  Introduction 

One  of  the  most  common  causes  of  significant  wind  shear 
is  the  gust  front  associated  with  thunderstorms.  The  thunder- 
storm gust  front  is  believed  responsible  for  the  Eastern  66 
crash  in  New  York  in  1975,  and  for  the  BOAC  crash  in  Kano, 
Nigeria  in  1965,  among  others.  The  severe  wind  shear  accom- 
panying thunderstorms  is  generated  by  a vigorous  rain-cooled 
downdraft,  which  spreads  out  horizontally  from  the  storm  cell 
as  it  approaches  the  ground.  The  cold  outflow  is  led  by  a 
strong  gusty  wind  which  often  occurs  as  much  as  20  km  ahead 
of  the  storm,  called  the  gust  front.  During  the  passing  of 
the  front,  winds  may  increase  from  a relative  calm  to  signi- 
ficant values  in  a very  short  time  and  then  decrease  as 
suddenly.  This  sudden  increase  in  wind  is  called  the  gust  size 
and  the  time  over  which  it  occurs  is  called  the  gust  length. 

4.2  Storm  Cells 

The  salient  kinematic  features  of  a thunderstorm  are  best 
described  by  considering  a particular  storm  reported  by  Kropfli 
and  Miller  [4-1],  The  storm  extends  to  a height  of  approxi- 
mately  13  km  and  covers  an  area  of  approximately  324  km  . 

Figure  4-1  schematically  illustrates  four  typical  features 
common  to  most  thunderstorms  which  are:  mid-level  (6  km  in 

i* 

this  case) , entrainment  of  cool  dry  air  from  behind  the 
storm,  a precipitation-filled  downdraft,  a gust  front 
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consisting  of  air  originating  from  the  middle  level  and  a warm 
air  updraft  ahead  of  the  storm.  Figure  4-2  from  [4-1]  shows 
the  approximate  three-dimensional  motion  of  air  parcels  in 
the  storm.  This  figure  illustrates  that  the  gust  front  at 
the  surface  on  the  right  forward  flank  has  its  source  of  air 
from  the  loft  rear  flank  and  exits  at  the  surface  under  the 
updraft.  The  source  of  air  for  the  northern  portion  of  the 
gust  front  comes  mainly  from  air  entrained  at  the  left  for- 
ward flank  and  partially  from  the  left  rearward  flank.  Air 
entering  at  mid-level  on  the  right  rear  flank  descends  and 
exits  at  the  rear  about  0.5  km  above  the  ground.  This 
downdraft  cut  the  backside  is  also  supplemented  with  inflowing 
air  which  starts  at  the  lower  edge  of  the  updraft.  Air 
entering  from  the  right  side  and  the  right  forward  flank  moves 
upward  with  cyclonic  curvature  and  exits  near  the  11  km  level. 
Most  of  the  upward  moving  warm  air  inflow  from  the  forward 
central  edge  of  the  storm  enters  at  a 25°  backward  tilt  (the 
direction  of  this  sheet  of  air  is  often  called  the  gust 
front  slope)  and  ascends  with  little  rotation,  exiting  near 
the  11  km  level.  Part  of  this  air  splits  off,  however,  and 
exits  from  the  left  side  near  the  9 km  level. 

Contours  of  vertical  velocity  relative  to  the  storm  motion 
in  horizontal  planes  through  the  storm  at  three  elevations 
are  shown  in  Figure  4-3.  Three  kilometers  above  the  surface 
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RELATIVE  TO  THE  STORM  MOTI 


the  inflow  air  is  moving  upward  at  approximately  3 m/s 
whereas,  the  downdraft  descends  in  the  center  at  4 m/s. 

At  7 km  the  inflow  has  penetrated  the  storm  with  vertical 
velocities  reaching  11  m/s.  The  highest  updraft  reported 
during  the  life  of  the  storm  occured  near  the  top  and  was 
observed  to  have  attained  a speed  of  19  m/s.  The  nature  of 
the  relative  velocity  in  the  horizontal  planes  is  indicated 
in  the  hodograph  shown  in  Figure  4-4.  These  velocity  contours 
correspond  to  a particular  instant  of  time  and,  although  re- 
taining similar  features,  will  show  strengthening  or  weakening 
as  the  storm  matures  and  decays. 

The  foregoing  comments  were  intended  to  describe  the 
large  scale  characteristics  of  a thunderstorm.’ 

It  is  the  air  motion  in  the  gust  front  of  the  thunder- 
storm below  the  2 km  level  which  is  of  most  significance  to 
aeronautical  operation  in  terminal  areas.  Some  relatively 
detailed  wind  tower  measurements  cf  gust  fronts  are  reported 
[4-2,  3,  4,  5]  although  a complete  understanding  of  the 
physical  mechanism  and  a predictive  mathematical  model  of  the 
three-dimensional  motion  is  still  needed. 

4.3  Gust  Fronts 


4.3.1.  Gust  Front  Types 

Figure  4-5  illustrates  the  common  features  of  a gust 

front. 

The  figure  shows  that  the  pre-storm  warm  moist  air  is 
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FIGURE  4-5  COMMON  FEATURES  of  a gust  front 


displaced  upward  by  the  more  dense  outflow  from  the  thunder- 
storm. The  mechanical  forcing  occurs  in  a 1 km  wide  band 
in  advance  of  the  gust  front.  After  rising  about  800  m 
over  the  head,  gravitational  instability  forces  the  uplifted 
air  down  in  the  wake  of  the  gravity  current  and  horizontal 
divergence  assists  entrainment  of  warm  air  across  the  in- 
version boundary  into  the  cold  outflow.  This  wake,  then, 
is  a highly  turbulent  zone  frequently  characterized  by  large 
shears  in  the  horizontal  wind  and  large  oscillations  in 
vertical  motion.  Similar  characteristics  are  observed  in 
secondary  outflow  surges.  Horizontal  divergence  has  been 
observed  to  be  so  strong  in  the  wake  of  these  secondary 
surges  that  large  volumes  of  air  may  descend  almost  to  the 
surface . 

The  extent  and  intensity  to  which  the  features  illustrated 
in  Figure  4-5  may  be  present  depend  upon  the  evolution  c.i  the 
thunderstorm.  Goff  [4-2]  has  grouped  his  data  into  four 
frontal  cases.  These  are: 

1.  Gust  fronts  associated  with  intensifying  storms  or 
accelerating  outflow. 

2.  Gust  fronts  associated  with  mature  intense  storms 
or  strong  outflow. 

3.  Gust  fronts  associated  with  dissipating  storms  or 
outflow  decelerating  with  respect  to  the  storm. 

4.  Gust  fronts  in  the  final  stage  of  life  cycle. 
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4.3.2  Gust  Front  Speed 

Table  4-1  lists  the  gust  front  speed  of  the  17  thunder- 
storms reported  in  [4-2’ . Gust  fronts  associated  with 
intensifying  and  accelerating  outflow  tend  to  move  more 
slowly  relative  to  the  ground  than  those  associated  with 
storms  in  the  state  of  final  dissipation  which  move  more 
rapidly.  The  gust  front  speed  has  been  related  to  the 
velocity  component  in  the  cold  air  by  Clarke  [4-61  as: 

c *=  0.67  uc  4.1 

Frank  [4-5]  gives  c = 0.62  uc  and  Cclmer's  [4-3]  data  from 
a visual  inspection  seem  to  support  a value  on  the  order  of 
0.67.  Goff  [4-2]  finds  Equation  4.1  suitable  for  case  1 
type  storms , but  that 

c = 0.7  uc  + 0.3  uw  4.2 

v?here  uw  is  the  velocity  component  in  the  warm  air  (normally 
negative)  correlates  types  2 and  3 storms  somewhat  better. 
.Type  4 storms  appear  to  have  a different  propagation 
mechanism  and  are  not  well  predicted  by  either  Equation  4.1 

or  4.2. 

4.3.3  Gust  Front  Characteristics 

Most  all  gust  front  passages  are  either  dry,  or  almost 
coincident  with  the  onset  of  rain.  The  dry  gust  front  may 
move  as  far  as  10  to  12  km  ahead  of  the  precipitation  and 
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TABLE  4-1 


GUST  FRONT  SPEEDS  OF  17  DIFFERENT  STORMS 


Type 


No.  of 
Storms 


Gust  Front 
Speed  (m/s-1) 


6.1,  5.0,  8.6,  11.6 


16.7,  11.0,  11. B,  8.5 
11.5,  13.1,  12.4,  5.5 


9.6,  8.0,  11.4,  6.9 


17.5 


can  travel  this  distance  from  the  cell  center  in  15  to  20  min, 


[4-7] . Goff  [4-2]  reports  that  the  separation  distance 


between  precipitation  and  the  gust  front  is  largest  for 


mature  storms,  case  2,  whereas  Colmer  [4-3]  reports  that  for 


the  11  thunderstorms  he  studied,  the  gust  fronts  remain  be- 


teen  0 to  7 km  from  the  storm  center  while  the  storm  is  in 


the  mature  and  developing  phases  and  then  move  ahead  of  the 


storm  by  8 to  20  km  during  the  decay  phase.  Goff  [4-2]  re- 


ports that  the  storm  overruns  the  gust  front  in  the  decaying 


stages  and  precipitation  tends  to  coincide  with  the  gust 


front  passage.  It  appears  that  the  definitions  of  the  storm 


life  by  these  two  authors  may  be  inconsistent  or  that  the 


storms  are  different.  Note  the  average  separation  distance 


for  case  1 storms,  from  Goff,  is  7.6  km;  case  2 storms  is 


10  km  which  correspond  roughly  to  Colmer *s  0-7  km  and  8-20  km 


range , 


^ ^ Vi- 


The  important  observation  from  the  above  is  that  gust 
fronts  with  sudden  wind  changes  may  pass  over  an  airport 
several  miles  ahead  of  any  indicated  precipitation  or  storm 
activity.  This  creates  unexpected  wind  hazards. 

4.3.4  Gust  Front  Shape 

The  gust  front  boundary  is  defined  as  the  boundary 
separating  the  flow  originating  from  a cold  source  and  the 
flow  of  air  originating  from  a warm  air  source.  The  boundary 
is  generally  shaped  with  a nose  protruding  ahead  of  the  front. 

Surface  drag  tends  to  retard  the  outward  rushing  cold  air 
near  the  ground  causing  .it  to  overrun  itself  above.  The 
point  of  the  nose  is  on  the  order  of  100  to  300  m above  the 
surface  (both  Colmer  [4-3]  and  Frank  [4-5]  report  the  peak 
at  250  m) , and  can  extend  horizontally  as  far  as  200  m into 
the  warm  air. 

The  lag  of  ground  level  air  behind  the  foremost  part 
of  the  front  is  expected  to  be  intermittent  with  the  cold 
air  falling  through  the  warm  air  entrained  beneath  the  nose 
and  again  reforming  due  to  surface  stress.  This  over- 
turning implies  high  turbulence  in  the  nose  region.  The 
bulging  gust  front  is  most  prevalent  in  mature  storms. 

Above  the  peak  of  the  nose  the  gust  front  boundary 

slopes  backward.  Goff  gives  angles  of  45°  for  case  1 storms 

and  45°  to  75°  for  case  2 storms.  A time-height  plot  of 

streamlines  shown  in  Figure  4-6  clearly  illustrates  this 
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slope  which  coincides  with  the  zero  streamline.  As  previously 
mentioned,  Kropfli  and  Miller  [4-1]  report  a 25®  slope  which 
corresponds  to  the  decaying  stages  of  their  storm.  Colmer 
[ 4-3  ] notes  that  the  frontal  slope  is  not  linear  and  parts 
of  the  gust  front  will  have  much  larger  slopes  than  the 
average  (see  Figure  4-7).  He  correlates  the  ratio  of  maximum 
slope  over  a 30  sec.  period  to  the  average  slope,  R,  with 
the  gust  length,  H,  by  the  relationship 

log  R = 0.667  log  H + const,  4.3 

Depending  on  the  relative  position  of  the  storm  to  the  air- 
port, a landing  aircraft  will  encounter  the  most  severe  wind 
shear  at  different  heights  above  the  ground.  Also,  the  rain 
typically  follows  some  distance  behind  the  gust  front  and 
the  pilot  may  or  may  not  land  in  rain.  Since  for  simulation 
purposes  it  is  useful  to  relate  the  location  of  the  shear 
and  the  probability  of  rain  to  the  position  of  the  airport, 
Figure  4-8  has  been  prepared  from  limited  data. 

Figure  4- 8 illustrates  the  altitude  at  which  an  air- 
craft approaching  an  airport  on  a 3°  glide  slope  might 
expect  to  encounter  a gust  front  which  has  passed  over  the 
airport  at  the  time  indicated  on  the  horizontal  scale.  The 
probability  of  rain  at  the  airport  relative  to  the  position 
of  encountering  shear  is  also  shown.  As  an  example,  the 
curve  shows  that  if  the  gust  front  had  passed  the  airport 
10  min.  prior  to  the  airplane's  approach,  the  pile*- 
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FIGURE  4-7  NONLINEAR  FRONTAL  SLOPE  AS  DESCRIBED  BY  COLMER,  [4-3]. 


ALTITUDE  AT  WHICH  AN  AIRCRAFT  APPROACHING  AN 
AIRPORT  MIGHT  EXPECT  TO  ENCOUNTER  A GUST  FRONT 
WHICH  PASSED  THE  AIRPORT  EARLIER. 


expect  shear  at  210  m and  a 33%  chance  of  landing  in  rain. 
This  curve  is  interpolated  from  the  data  of  the  preceeding 
discussion  and  should  be  used  only  as  a rough  estimate. 

4.3.5  Gust  Characteristics 

The  gust  size  (magnitude  of  wind  increase)  and  the  gust 
length  (length  of  time  over  which  the  wind  increases)  were 
plotted  by  Colmer  for  11  case  studies.  On  an  average,  the 
gust  size  increases  by  50%  between  the  surface  and  500  m 
with  the  most  increase  occurring  in  the  first  50  m.  This  is 

a significant  observation  since  a ground  based  anemometer 
will  read  winds  considerably  less  than  those  a hundred  or 
so  meters  above.  ? 

Note,  however,  that  the  gust  size  profile  exhibits 
considerable  variability.  The  gust  length  on  the  other 
hand  shows  less  variability  and  an  almost  linear  decrease 
with  height.  Additional  gust  size  and  length  data  are 
given  by  Sinclair  et  al  [4-4]  for  three  Florida  storms. 

They  report  at  18  m an  average  gust  size  of  4.8  m/s  and  a 

gust  length  of  12.8  s,  and  at  the  height  where  the  peak  wind 

occurs,  they  report  an  average  gust  size  of  6.6  m/s  and  a 

gust  length  of  14.8  s.  The  length  scale  based  on  a constant 

gust  front  speed  (which  was  not  reported)  does  not  show  the 
same  decrease  with  altitude  that  Colmer  reports.  Goff  [4-2] 
finds  the  change  in  the  horizontal  wind  component  normal  to 
the  gust  front  to  have  a value  between  4 and  20  m/s  with  an 
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average  of  12.8  m/s  which  shows  good  agreement  with  Colmer's 
average  value  above  50  m. 

The  track  of  the  maximum  gust  tends  to  coincide  with  the 
track  of  the  storm  cell.  The  size  of  the  gust  decreases 
almost  linearly  from  this  maximum  value  havinc  an  approximate 
15%  decrease  in  gust  size  5 km  from  the  center  and  a 32% 
decrease  10  km  from  the  center. 

Frequently  more  than  one  gust  is  experienced  during  the 
passage  of  a thunderstorm.  Ward  and  Arnett  [4-8]  report 
multiple  cold  air  surges  for  a single  thunderstorm  which 
they  attribute  to  storm  pulsation.  Goff's  [2-2]  data  also 
shows  secondary  surges.  There  appears  to  be  no  reduction 
of  the  data  which  allows  a quantitative  statement  about 
these  multiple  surges  to  be  made.  Goff  [4-9]  in  a recent 
paper  discusses  secondary  surges  of  high  momentum  air  which 
appeared  in  the  outflow  air  mass  from  the  thunderstorm  on 
June  7,  1971.  Two  surges  observed  were  both  characterized 
by  a strong  shear  zone  in  the  Wx-component  and  large  vari- 
ations in  the  vertical  wind  speed.  In  one  surge,  a downdraft 
in  excess  of  11  m/s  was  observed  at  the  177  m level. 


4.3.6  Wind  Speed  Fields 

The  most  useful  data  for  simulation  of  aircraft  flights 
through  gust  fronts  is  the  three-dimensional  velocity  field 
in  the  vicinity  of  the  gust  front  as  a function  of  position 
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and  time.  A step  toward  providing  these  data  is  the  height- 
time plots  given  by  [4-2,  4,  5].  Often  these  plots  are 
converted  to  height-distance  plots  by  assuming  Taylor's 
hypothesis.  Taylor's  hypothesis  allows  wind  speeds  measured 
in  time  at  a point  in  space  to  be  converted  to  wind  speed 
variation  with  position  along  a line  in  the  direction  of  the 
mean  wind  by  assuming  the  turbulent  flow  is  carried  past  the 
point  of  measurement  sufficiently  rapid  that  the  turbulence 
pattern  is  "frozen"  in  the  air  flow  and  does  not  change  over 
the  time  the  measurement  is  taken  (ie  position  = mean  velocity 
X time) . Typical  height-time  plots  are  shown  in  Figure  4-9 
and  height-distance  plots  in  Figure  4-10.  Generally,  these 
data  are  highly  smoothed  by  filtering  techniques  which  remove 
the  high  frequency  content.  Whereas,  Figure  4-9  gives  the 
wind  speed  and  the  wind  direction  contours  in  the  vertical 
plane.  Figure  4-10  expresses  the  data  in  terms  of  the  three 
velocity  components:  W , vertical  velocity,  W , horizon.  ..1 

velocity  normal  to  the  front,  and  W , lateral  wind  speed 
parallel  to  the  front. 

The  data  corresponds  with  the  schematic  of  the  gust 
front  shown  in  Figure  4-3.  Inspecting  the  component  of 
wind  speed  in  Figure  4-10  one  observes  a negative  wind  speed 
corresponding  to  warm  air  inflow  toward  the  front  which  is 
delineated  by  the  zero  wind  speed.  The  positive  values  in 
the  height-time  plots  correspond  to  cold  air  from  the  thunder- 
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storm  flowing  toward  the  front.  Negative  wind  speed  near  the 
surface  in  the  cold  air  zone  indicates  an  undercurrent  prob- 
ably resulting  from  the  cold  air  riding  over  a low  level 
temperature  inversion.  Although  filtered  out  of  the  data 
of  Figure  4-10,  a highly  turbulent  wake  region  expanding  in 
depth  with  increasing  distance  behind  the  front  occurs  behind 
the  elevated  "head"  of  dense  air  which  is  characterized  by 
strong  turbulence  and  mixing. 

Inspection  of  the  vertical  wind  speed  in  Figure  4-10 
shows  a vertical  updraft  attributed  to  lifting  of  the  less 
dense,  warm  air  over  the  more  dense  forward-moving  air  mass. 
Maximum  updrafts  generally  occur  1 km  or  less  ahead  of  the 
gust  front.  The  magnitude  of  these  updrafts  generally  in- 
creases wi*-h  frontal  strength  (increasing  frontal  slope)  . 

The  value  of  the  maximum  updraft  reported  by  Goff  [4-2]  is 

6.7  m/s,  and  that  by  Browning  and  Harold  [4-10]  is  In  excess  of 
8 m/s  and  10  m/s  for  two  different  fronts.  The  record  from 
the  digital  flight  data  recorder  of  Eastern  Air  Lines  Flight 

66  and  Eastern  Air  Lines  902  for  the  thunderstorm  over  J.F. 
Kennedy  Airport  gave  updrafts  of  approximately  4.3  m/s,  and 
5.2  m/s,  respectively.  However,  these  values  probably  do  not 
correspond  to  the  maximum. 

4.3.7  Wind  Shear 

Strong  vertical  and  horizontal  shears  occur  with  all 
gust  fronts.  Vertical  shear  is  a variation  in  the  wind 
speed  components  with  height,  z,  whereas  horizontal  shear 
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refers  to  variation  with  horizontal  distance  (in  this  case, 
WAtli  the  longitudinal  mean  wind  direction,  x)  . Large 
vertical  shear  (3()/3z)  occurs  near  the  surface  and  at  the 
upper  boundary  and  large  horizontal  shear  (D()/3x)  occurs 
at  the  outflows  leading  edge.  Secondary  horizontal  shear 
may  occur  also  due  to  multiple  outflow  surges. 

The  maximum  wind  shear  averaged  from  the  20  thunderstorm 
cases  studied  by  Goff  f 4—2 ] are  tabulated  in  Table  4-2a  and 
Table  4-2b.  The  spatial  derivatives  were  formed  from  the 
relationships  3/3x  = c~^3/3t.  From  the  table  one  notes  that 
shears  of  horizontal  wind  are  greatest  near  the  ground 
reaching  a minimum  at  200  m and  then  increasing  upward.  This 
tends  to  conflict  with  Colmer's  [4-33  conclusion  that  due  to 
the  increase  of  gust  size  and  decrease  of  gust  length  with 
altitude,  horizontal  wind  shear  in  the  gust  front  between 

the  surface  and  500  m above  ground  doubles,  i.e.,  the  average 
(of  11  cases)  surface  shear  is  0.0063  s ^ while  the  500  m 

level  shear  is  0.016  s It  should  be  noted  that  reported 
wind  shear  values  are  based  on  a very  sparse  grid  system  and 
on  highly  smoothed  data  and  are,  therefore,  at  best  gross 
approximations . 


4 . 4 Wind  Shear  Prediction  Techniques 


4.4.1  Raw  Data 

Despite  the  filtering  and  course  grid  space  used  to 
obtain  the  spatial  distribution  of  wind  speeds  and  wind 
shears  given  in  the  previous  section,  the  results  are  ex 
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TABLE  4- 2 a 

Maximum  Absolute  Shear  (m  s~^  100  m-^)  [4-2] 

*Grid  level  or  layer  of  gradient: 

1 2 3 4 5 6 7 8 9 10 


8Wx/  8x 

11.1 

5.0 

4.9 

5.4 

4.7 

4.1 

4.8 

5.3 

4.5 

4.6 

aw/  ax 

7.0 

3.2 

3.4 

3.8 

3.4 

3.7 

4.0 

4.8 

4.8 

7.8 

awz / ax 

0.0 

1.2 

2.2 

3.4 

3.8 

3.4 

2.6 

2.8 

3.3 

3.8 

awx/  3z 

16.2 

8.4 

8.7 

9.6 

9.4 

8.6 

10.8 

6.5 

10.6 

— 

awy/  3z 

26.1 

10.2 

6.9 

5.6 

8.2 

6.4 

6.4 

6.3 

9.7 

— 

3WZ/  3z 

3.8 

2.1 

2.3 

2.0 

2.6 

3.1 

3.1 

3.1 

3.1 

— 

TABLE  4 -2b 

Average  Maximun 

Absolute  ! 

Shear 

(m  s"1 

100  m" 

1)  [4-2] 

* Grid  level  i 

or  layer  of  gradient: 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

3W  / 3x 

X 

4.3 

2.0 

2.3 

2.3 

2.3 

2.4 

2.6 

2.9 

2.9 

3.2 

3W  / 3x 

V 

4.0 

1.5 

1.5 

i . 7 

1.9 

1.9 

1.9 

2.2 

2.2 

2.6 

3Wz/3x 

0.0 

0.6 

0.9 

1.2 

1.3 

1.2 

1.1 

1.5 

1.9 

2.4 

3W  / 3z 

X 

11.6 

5.8 

4.7 

4.6 

4.1 

4.7 

5.0 

4.1 

5.1 

— 

3W  / 3z 
y 

10.6 

5.3 

4.1 

3.6 

3.9 

3.6 

3.4 

3.2 

4.3 

— 

3Wz/3z 

2.1 

1.6 

1.7 

1.3 

1.7 

1.9 

1.9 

1.9 

1.9 

— 

* Levels  1 through  10  are  50  meters  apart. 
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tremely  useful  for  estimating  the  influence  of  gust  fronts 
on  aircraft  flights.  Two  ways  of  using  the  data  are  to 
store  them  in  the  computer  on  an  appropriate  grid  system 
and  use  linear  interpolation  with  a table  look-up  scheme 
to  give  the  wind  input  values  to  the  computerized  airplane 
equations  of  motion;  or  to  superimpose  geometrically  scaled 
flight  paths  on  to  the  contour  plots  and  pick  off  the  winds 
"seen"  by  an  aircraft  traversing  that  flight  path.  The 
latter  method  is  more  restrictive  in  that  the  aerodynamic 
forcing  functions  computed  from  the  input  winds  assume  the 
airplane  remains  on  or  in  close  proximity  to  the  selected 
flight  path. 


4.4.2  Mathematical  Models 

Mathematical  schemes  for  computing  wind  fields  associ- 
ated with  thunderstorm  gust  fronts  are  still  in  the  formative 
stages.  The  most  direct  mathematical  approach  is  probably 
that  of  Mitchell  [4-11]  who  solved  the  finite  difference 
equation  for  the  two-dimensional,  time  dependent  primitive 
equation  of  fluid  motion.  Surface  friction  is  modeled  with 

a bulk  aerodynamic  drag  function  applied  at  the  lower  levels 
of  the  computation  grid.  Vertical  diffusion  of  momentum  is 
modelled  with  linear  diffusion  terms,  and  velocity-slip  is  ' 
permitted  at  the  surface. 

This  model  appears  to  capture  much  of  the  essential 
physics  of  the  flow  and  provides  a computational  tecnnique 
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which  might  be  incorporated  into  a computer  code  for  simul- 
taneous solution  of  the  fluid  mechanics  with  the  airplane 
dynamics.  Unfortunately,  however,  the  computer  solution 
is  costly  and  this  prohibits  any  extensive  flight 
simulation. 

A less  sophisticated  approach  to  modelling  two-dimensional 
flow  fields  about  a gust  front  is  proposed  by  Sinclair  et  al 
[4-4],  They  scaled  the  three  Florida  thunderstorms  with  the 
scaling  parameters;  time,  At,  which  characterizes  gust  length; 
velocity  AV,  which  characterizes  gust  size  and;  height,  Az, 
which  is  the  height  of  maximum  velocity.  Figure  4-11  shows 
scaled  velocity  contours  which  Sinclair  et  al  [4-4]  conclude 
arc  quite  similar  to  all  other  scaled  gust  profiles  and, 
hence,  if  AV,  At,  and  Az  are  known,  at  least  part  of  the 
structure  of  the  wind  front  can  be  determined.  Because  of 
the  difficulty  in  evaluating  Az,  this  approach  was  abandoned 
and  a simple  hypothesis  to  find  the  wind  structure  was 
adopted. 

We  will  not  pursue  this  approach  further  since  the 

effort  of  reference  [4-4]  was  toward  forecasting  gust  fronts, 
whereas  our  goal  is  to  provide  realistic  wind  models  for  use 
in  parametric  simulation  and  design  studies.  The  scaled  gust 
fronts  of  Figure  4-11  provides  such  a model  in  that  it  allows 
us  to  conduct  design  analysis  and  simulation  studies  with 
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parametric  variation  of  gust  length,  gust  size  and  height 
of  maximum  velocity.  Since  we  wish  to  establish  design 
envelopes  or  limits  rather  than  forecast  gust  fronts,  we  are 
at  liberty  to  select  any  values  of  Az,  AV,  and  At  within  the 
confines  of  good  engineering  judgement.  More  will  be  described 
of  this  approach  in  the  final  report.. 

Updrafts  and  downdrafts  associated  with  gust  fronts  along 
a given  approach  path  have  been  correlated  by  Fichtl  and 
Camp  [4-10].  They  scaled  vertical  wind  speeds  along  a 3° 
glide  slope  from  the  gust  front  data  of  Goff  [4-2]  and  also 
incorporated  the  vertical  wind  speeds  reconstructed  from  the 
digital  flight  data  record  of  Eastern  66.  Their  mathematical 
representation  of  the  sequence  of  vertical  wind  speeds  en- 
countered by  an  aircraft  during  landing  is  given  by  the 
fol lowing : 

Major  Downdraft: 
x-x 


w 


-P]A  sin|»-qy ].  ~D  ... 


>z  >z 


4.4 


Major  Updraft: 

(l-2qo)  (x-xr)3  + (l"3q^)  (x-xr)2  + [2q^3^  . 


w = A 


z > z > z - L 
r - r 


4 . 5 
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Minor  Downdraft: 

x -1  -x 

w - -p2A  sin  [it  — ] ; zr  - L > z ^ zr  - (1  + q2)L  4.6 


Minor  Updraft: 

x-l-x 

w * +P2A  sinliT  

where 


; zr-  (1  + q2)  L>  z>  zr  "(1  + 2q2)L 

4.7 


z 

L 


4.8 


The  various  quantities  in  the  above  equations  are  defined 
as  following: 

w = thunderstorm  cold  air  outflow  vertical  velocity 
z = altitude  of  airplane  above  surface 
Zp=  altitude  of  the  top  of  the  major  downdraft 
zr=  altitude  of  the  top  of  the  major  updraft 
Zg*=  altitude  of  the  bottom  of  the  minor  downdraft 
zm=  altitude  of  maximum  updraft  velocity  in  the  major 
updraft 

A = amplitude  of  major  vertical  velocity  updraft 
L = vertical  extent  of  major  vertical  velocity  updraft 
relative  to  the  flight  path 

p^=  ratio  of  major  downdraft  to  major  updraft  velocities 
P2=  ratio  of  minor  downdraft  or  minor  updraft  to 
major  updraft  velocity 

V (Vzm)/L 


*This  parameter  is  not  required  in  the  calculation 
but  serves  to  physically  define  qQ. 


9l  * (2r'zDIA 
<32  = (zr‘“zB"’k) /k 

Values  for  the  cold  air  outflow  parameters  are  given  in 
Tables  4-3  and  4-4.  A schematic  is  given  in  Figure  4-12. 
Table  4-3  provides  cold  air  outflow  parameters  for  typical 
vertical  wind  speeds  as  derived  from  data  provided  by  the 
NOAA/NSSL.  Table  4-4  provides  cold  air  outflow  parameters 
which  encompass  the  worst  case  (Eastern  66) . 

To  apply  the  model,  one  enters  the  model  with  the 
parameters : 

zr'  zd'  V ^1 ' ^2 

consistent  with  the  range  of  parameters  in  Tables  4-3  and 
4-4.  It  is  recommended  that  a family  of  profiles  which 
encompass  the  full  range  of  variation  of  the  model  as  implied 
in  Tables  4-3  and  4-4  be  used.  This  will  insure  that  certain 
conditions  will  not  be  overlooked. 

Wind  speed  as  a function  of  height  from  combined  tower 
and  rawinsonde  data  near  the  time  of  the  wind  speed  maximum 
was  correlated  by  Frank  [4-5]  with  the  expression 

W = wn+  tan  ( (z-z  )/  AH)  (l-e"z/zm)  4.9 

X v o 
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where 


V 

94. 

72m 

AH= 

889 

m 

V 

402 

m 

o* 

n 

25. 

0 m/s 

AW= 

-23 

. 5 m/s 

Figure  4-13  shows  the  results.  Also  Frank  [4-5]  notes  that 
in  the  surface  boundary  layer  regions  the  wind  profile  obeys 
a power  law 


xr 


(z/zr)n 


4.10 


where  W is  the  wind  speed  normal  to  the  front  at  height  z, 
zr  is  a reference  level,  and  Wxr  = Wx(zr).  Just  behind  the 
front  n was  found  equal  to  0.39  from  tower  data  and  equal 
to  0.35  from  Doppler  radar  data.  The  power  law  was  reported 
to  give  gpod  agreement  with  gust  front  data  up  to  350  m. 

Sinclair  et  al  [4-4]  demonstrated  that  a logarithmic 
wind  law  is  valid  from  the  surface  to  approximately  100  m. 
Therefore,  the  ratio  of  wind  speed  at  two  different  levels 
becomes  simply  a function  of  surface  roughness,  i.e.. 


V2  ln(z2/zQ) 

“ InTYjTz^) 


4.11 
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TABLE  4-3 

Typical  Values  of  Cold  Air  Outflow  Parameters  (Based  on  NOAA/NSSL 
Data)  [4-10] 


Parameters 


Units 


64  + 0.223zr  <L  <183  + zr 
0<  A < 1 + 0 . 011z„ 


ft  sec 


0 . 3<  pL  < 1.2 
0<p2£  0. 35 
qQ  = 0.36 

0.25  qL<  1 

°1  ^2—  2-3 

L<  z < 1500 

- r— 

ZD  = Zr  " <*1L 
zB  = zr-  (1  + q2)L 


Nondimensional 

Nondimensional 


Nondimensional 


Nondimensional 


Nondimensional 


TABLE  4-4 


Worst  Case  Values  of  Cold  Air  Outflow  Parameters  (Based  on 
Eastern  66  and  NOAA/NSSL  Data)  [4-11] 


Parameters 

64  + 0.223  zr  < L < 183  + zr 

0 < A £ 10  + O.Ollz 

0-3  < px  < 1.2 

0 1 p2  < 0. 35 

q0  = 0.36 

0 < q 1 < 2 

05  q2  1 2.3 


L < zr  < 1500 
zD  = zr  - q:L 


Units 


ft  sec 


Nondimensional 

Nondimensional 

Nondimensional 

Nondimensional 

Nondimensional 


Zg  zr~  (1  + q2)  L 
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FIGURE  4-12  VERTICAL  WIND  SPEED  AL0N6  A 3*  GLIDE  SLOPE 
THROUGH  GUST  FRONTS  AS  CORRELATED  BY 
F1CHTL  AND  CAMP,  [4-10]. 
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5.0  BOUNDARY  LAYER  OVER  FLAT  TERRAIN 

5.1  Introduction 

Considerable  information  on  the  behavior  of  the  atmo- 
spheric boundary  layer  at  specific  locations  has  been  obtained 
with  individual  meteorological  towers  up  to  altitudes  of 
typically  150  m with  some  data  up  to  500  m.  These  data  are 
normally  assumed  valid  over  flat,  horizontally  homogeneous 
terrain  and  are  transferred  to  spatial  coordinates  with 
Taylor's  hypothesis  which  requires  that  the  boundary  layer 
be  satistically  stationary.  Above  the  altitude  which  can  be 
reached  with  towers;  aircraft,  rawinsonde  and  jimsphere  data 
provide  insight  into  the  turning  layers  of  the-  atmosphere 
but  these  data  are  not  as  precise  as  tower  data. 

5 . 2 Mean  Flow  in  the  Lower  Surface  Layers 

Below  150  m the  mean  wind  direction  has,  in  general, 
little  variation  with  height,  and  over  flat  terrain  is  con- 
sidered to  be  steady  plane-parallel  flow  of  a horizontally 
homogeneous  and  statistically  stationary  surface  layer. 
Numerous  reviews  of  the  atmospheric  boundary  layer  are 
available  [5-1,  2,  3,  4,  5]  and  it  is  only  necessary  to 
summarize  these  here. 

In  the  surface  boundary  layer,  the  mean  wind  is  con- 
sidered to  be  influenced  by  surface  conditions,  atmospheric 
stability,  and  elevations.  Mean  wind  models  receiving  the 
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greatest  acceptance  are  expressed  in  terms  of  nondimensicmal 
wind  shear  as: 

zk  dW 

— 37"  “ *(Z/L) 

The  wind  speed  is  then  given  by 

z-fz 

V>  J tifZU  d2 

z 
o 

W (z)  * 0 

Wz(z)  - 0 

where  Wx(z)  , V.Ty(z)  and  W2(z)  are  the  wind  speed  in  the  x, 
y,  and  z directions  respectively;  the  x axis  is  aligned 
with  the  mean  wind  direction,  the  y axis  is  perpendicular 
to  x and  the  z axis  is  vertical  to  the  horizontal  plane. 

The  constant  k is  the  von  Karman  constant  having  a numerical 
value  of  approximately  0.4  and  zQ  is  a length  scale  that 
characterizes  the  surface  roughness.  Typical  values  of  zQ 
are  given  in  Figure  5-1. 

The  surface  friction  velocity  u*  is  normally  assumed 
constant  which  is  experimentally  verified  only  in  the  lower 
30  meters,  however,  its  variations  to  150  m is  generally 
negligible . 
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5.2 

5.3 

5.4 
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Richardson's  number  is  a nondimensional  ratio  between  the 
mechanical  wind  shear  that  tends  to  displace  air  and  the 
buoyancy  force,  which  may  damp  or  amplify  this  tendency. 

The  atmosphere  is  said  to  be  stable  if  Ri,  z/L'>0  corre- 
sponding to  3T/3z<-g/Cp  indicating  a stable  lapse  rate  or 
inversion;  neutral  if  Ri,  z/L'=0  giving  3T/3z = -g/Cp 
indicating  an  adiabatic  lapse  rate;  and  unstable  if  Ri, 
z/L' <0  giving  3T/3z>-g/Cp  indicating  an  unstable  lapse  rate. 

The  empirically  derived  expressions  for  <f>(z/L')  for 
different  conditions  of  stability  are  given  by  a number  of 
investigators.  Barr  et  al  [5-5]  recommends  the  following 
forms  for  <f>(z/L'). 

Neutral  case: 

<p  (z/L' ) = 1;  z/L'  = 0 5.10 

Stable  case; 

<j> (z/L'  ) = 1 +a'£,  j 1>  z/L'  > 0-  5.11 

Very  stable  case: 

<f>  (z/L'  ) - 1 + o'  ; z/L'  >1  5.12 

where 

- 1 + o'  Ri 
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Unstfcole  case 
4>  ( z/L' ) * 


(1  - Y’  Ri) 


¥; 


0 >z/L' 


5.13 


where 


z 

L' 


Ri 

(1+  y'  Ri^ 


The  nondimensional  wind  shear  has  the  form  shown  in 
Figure  5-2  where  a * 4.5  and  y * 4a*  18.  Details  of  the 
development  of  Equations  5.10  to  5.13  are  given  in  [5-5/6]. 
The  integrated  form  of  Equations  5.10  to  5.13  is 


2 «f  2 

fin  — --  + <Kz/L')] 
zo 


5.14 


where 


tKz/L') 


z/L' 

<b  (z/L 1 ) -1 

z JV 


d (z/L ' ) 


5.15 


Before  considering  the  result  of  the  integration,  it  is 
noted  that  the  above  expressions  for  nondimensional  shear 
become  increasingly  less  accurate  at  altitudes  above  z > 90  m. 
This  is  due  to  the  fact  that  their  development  assumes  a 
constant  shear  stress  which  does  not  hold  at  high  altitudes. 
Barr  et  al  [5-5]  propose  an  extension  of  the  model  to 
altitudes  on  the  order  of  350  m by  allowing  the  friction 
velocity  to  vary  linearly  as 


u 


u (l-z/6) 
•>V0 


58 


5.16 


where  6 is  the  atmospheric  boundary  layer  thickness  given  by 

6 = u*o/  10. It  5-17 

The  parameter  f is  the  Coriolis  parameter.  The  non dimensional 
wind  shear  is  then  given  by 


, 3W 
kz  x 


u*o  32 


= (1  - z/6)  4>(z/L') 


5.18 


Integration  of  the  shear  gives 


u . z + z 

W = -£■  [In  

x k z„ 


I + !» 


5.19 


The  wind  profile  for  the  given  stability  cases  then  becomes; 


Neutral  case: 


Uj.  z + 

W = [in  - | ] ; z/L'  = 0 

x k z _ o 


5.20 
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Stable  case: 


Vt?  nn  lli 


« £ 4 iv  * ^ 


y + <*'  J-,  (1-  a/26)  ]; 


0<z/L'<  1 


Unstable  case: 


Ujl  - Z "f*  Z 
Wx  ‘ TT  f ln  —I—2 


f + ♦,<!•  ■ })i 


where 


*,  <!•  f) 


-z/L' 


^(L.)  (l-|  ) + |^-f  / i/;(Od£ 


where 


£ * z/L'  and  z/L'  <0 


Numerical  integration  of  Equation  5.23  is  carried  out  in 
[5-5]. 


Very  stable  case: 


U*  z + z 

Vtt 

o 


+ cx'(l+  In  -,)  + £ [1  + 0». 
L 6 


2Tz/l’) 


z/L'  >1 


5.24 
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The  very  stable  wind  profile  is  applied  above  the  alti- 
tude Zj|  where  z^/L,=x  1. 

Equation  5.24  must  be  used  with  discretion  since  it  is 
based  primarily  on  conjecture.  Little  is  actually  known 
about  the  very  stable  wind  profile,  the  shape  of  which  is 
largely  determined  by  the  absence  of  turbulence  above  z/L*>l. 
Very  near  the  ground  where  z/L*  is  less  than  unity  a turbulent 
layer  can  exist.  As  height  increases  beyond  z/L'=l  a non- 
turbulent  layer  can  occur  which  is  decoupled  from  the  lower 
layer.  Change  in  wind  direction  in  excess  of  45°  is  not 
uncommon  between  these  uncoupled  layers.  Leurs  [5-4]  suggests 
a discontinuous  wind  profile  for  very  stable  conditions. 

Below  a given  level  z.  he  assumes  either  a calm  W (z)  =0 
or  a logarithmic  profile.  Above  z^  he  considers  various 
constant  values  of  W (z) . The  wind  speed  change  across 

X 

is  taken  as  discrete  which  is  probably  not  physically 
meaningful  but  represents  a limiting  case. 

The  application  of  the  mean  wind  profiles  is  most 
frequently  based  on  a reference  wind  speed  at  a specified 
height.  Present  efforts  to  standardize  tower  heights  at 
airports  to  20  ft  suggests  20  ft  might  serve  as  a meaningful 
reference  height.  Given  a reference  wind  and  height,  u*Q 
can  be  determined  from 


u 


*o 


Vzref> 


'ref  + 


+ *,<”LT 


ref  ref 


)- 


'ref 


In 
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In  view  of  the  fact  that  6 is  a function  of  u*Q  an  iteration 
procedure  is  required  if  the  z/S  function  proposed  by 
Barr  et  al  [5-5]  is  retained.  The  effects  of  altitude 
variation  of  shear  stress,  ie  u*-  u^  (l-z/6) , is  negligible 
for  unstable  conditions  but  can  be  appreciable  for  stable 
conditions.  The  Richardson  number  or  L'  must  also  be 
known  at  the  reference  height;  or  a parametric  study  of 
variation  with  Ri  may  be  carried  out  depending  on  the 
design  problem. 

The  roughness  length, zq, should  be  specified  for  each 
particular  airport  and  may  be  a function  of  wind  direction 
(see  [5-7]  for  an  example  of  a surface  roughness  wind 
direction  rose) , wind  speed  (for  instance  tall  grass  and 
water  have  changing  roughness  characteristics  with  wind  speed) , 
season,  etc. 

A comparison  of  the  various  equations  given  above 
(see  [ 5-5 } ) indicates  that  stable  conditions  result  in  more 
severe  and  higher  winds  above  a certain  level  for  the  same 
reference  wind  speed.  However,  since  this  wind  is  propagated 
from  the  geostrophic  wind,  it  might  be  expected  that  greater 
stability  at  the  reference  height  will  be  associated  with 
lower  mean  wind  speeds  at  the  same  height,  and  consequently 
the  differences  in  the  shears  between  stable  and  unstable 
conditions  may  not  be  so  great. 
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5.3  Turning  Layer 


5.3.1  Introduction 

Above  the  constant  stress  layer  the  presence  of  mean 
flow  synoptic  scale  horizontal  pressure  gradients  and  the 
action  of  the  Coriolis  force  tend  to  produce  vertical  vari- 
ation in  the  vertical  transport  of  horizontal  momentum  and 
clockwise  turning  of  the  mean  flow,  looking  toward  the 
earth,  in  the  northern  hemisphere  and  counterclockwise 
turning  in  the  southern  hemisphere.  As  indicated  earlier, 
however,  experimental  evidence  [5-8,9]  appears  to  show  that 
in  the  height  interval  30  m < z < 150  m the  vertical  variation 
of  the  wind  direction  is  approximately  2°  - 4°.  which  is 
negligible  for  most  design  analyses.  Additionally,  by  allow- 
ing u*  to  vary  linearly  with  altitude,  Barr  et  al  [5-5]  argues 
that  Equations  5.20  to  5.24  reliably  predict  W^fz)  to  an 
elevation  of  approximately  300  m.  This  may  be  true  under 
barotropic  conditions,  but  if  baroclinicity  exists,  the 
turning  of  the  wind  may  become  important  as  low  as  150  m and 
should  be  considered  in  analyses  of  aircraft  flight  char- 
acteristics and  control  system  design. 

Most  solutions  for  the  wind  profile  in  the  turning  layer 
treat  the  barotropic  case  which  assumes  that  the  large-scale 
synoptic  horizontal  pressure  gradients  which  are  related  to 
the  geostrophic  wind  W G by 

„ _ 1 3?  . w = 1 5.26 

wxG”  " pf  3y  ’ yC  pf  3x 


64 


do  not  vary  in  the  vertical  direction.  Surveys  of  such 
solutions  are  given  in  15-1,  3,  10,  11  and  12]  . Reference 

r 5—12 ] summarizes  as  well,  analytical  solutions  which 
discard  the  barotropic  assumption  and  consider  specified 
height  variation  of  the  geos  trophic  wind,  i.e.»  the  baroclinic 
case.  The  results  of  these  solutions  will  be  summarized 
below,  foi  further  details  the  reader  should  see  the  stated 
references . 

5.3.2  Barotropic  Turning  Layer 

According  to  similarity  theory,  if  a variable  wind  is 
appropriately  scaled  then  its  profile  is  given  by  an  empirically 
determined  universal  function.  For  flow  far  from  the  surface, 
similarity  gives : 


W - W 


ref  _ 


u 


= F 


*o 


C-  -) 

V 1/ 


5.27 


where  ^vef  is  the  appropriate  wind  scale  to  be  discussed  later 
and  F is  the  outer  profile  function. 

For  flow  in  the  surface  layer 


W 

u* 


-*■ 

l 


5.28 
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TABLE  5-1 


Values  of  A(h/L)  and  B (h/L)  Under  Neutral  Conditions  h/L  - 0 


Re ference 

A (0) 

B (0) 

Yamada  [5-15] 

1.855 

3.020 

Fichtl  [5-1] 

1.6 

4.7 

Tannekes  [5-10] 

2 

5 

Monin  [5-3] 

1 to  2 

4.2  to  4.7 

Clarke  and  Hess  [5-16] 

1.1  ± 0.5 

4.3  ± 0.7 

and  the  angle  between  the  direction  of  the  surface  wind  and 
that  of  the  reference  wind  vector  Wref  is  obtained  from 


u* 

sin  a "fcffw — 


f 

TfT 


5.33 


For  the  barotropic  atmosphere  the  appropriate  scale  for 
wind  is 


ref 


WG 


5.34 


where  the  geos trophic  wind,  WG,  is  constant. 

The  appropriate  height  scale  for  neutral  stability  is 
h = c^u*/(f j where  c^  = 0.3.  Departures  from  neutral 
stability  cause  some  controversy  over  the  appropriate  form 
of  h.  Yamada  [5-15]  found  the  relationship  between  h/L  and 
u*/|f|L  shown  in  Figure  5-3.  The  expression  h/L  = 0.3  u*/|f|L 
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holds  reasonably  wall  for  u*/Jf|L<2QQ  but  it  underestimates 
h/L  considerably  for  u*/jf|L>  200.  Other  workers  [5-17] 
find  the  height  of  the  inversion  layer  z ^ during  nighttime 
and  the  mixing  layer  during  daytime  to  be  the  correct  scaling 
height. 

The  foregoing  analysis  does  not  permit  computation  of 

the  wind  profile  variation  with  altitude.  It  does  express 

the  geostrophic  drag  coefficient  u*/|WQ|  and  the  cross  iso- 

baric  wind  angle  as  a function  of  the  surface  Rossby  number 

Ro  = |w  |/fz  for  given  values  of  h/L.  Note,  if  h can  be 
6 o 

written  u*  /f.  Equation  5.33  can  be  rewritten 


|Wf|  k 2 1 Wp | 2 % 

In  Ro  * B (h/L)  + In  ^ • +[■— -JS'-  A (h/L)  2 ] 


5.35 


u 


*0 


u 


*0 


A meaningful  wind  profile  variation  with  elevation  can 
not  be  expressed  in  a simple  mathematical  form.  Most  all 
solutions  to  the  wind  spiral  with  height  have  been  carried 
out  numerical!  . The  equations  solved  are 


1 3Tx  . 


P 32 


- -f<WxG-V2)) 


and 


- f<v  vz)) 


5.36 


5.37 
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where 


dWx  dWv 

pK  di~  : Ty=  pK  ^ 


dz 


5.38 


The  eddy  momentum  K must  be  modeled  to  close  the  set  of 
equations.  The  expression 


9W  2 

K = a2  [ ) + (-L)  y 


8Wv  2 
Tz 


5.39 


where 

k(z  + zQ) 

1 ' 17  T-u/^575 


5.40 


was  proposed  by  Lettau  as  reported  by  Hanna  [5-11]. 

Other  authors  [5-18,  19]  have  used  similar  forms  for  K 
and  Z and  obtain  essentially  che  same  solutions  for 
W^(z)  and  W^(z).  Figure  5-4  shows  Lettau’ s solution  in 
dimensionless  form.  This  solution  is  valid  for  all  neutral 
wind  spirals, the  only  variation  being  the  length  of  the 
straight  line  extension  along  the  left  asymptote. 

To  construct  the  neutral  oarotropic  wind  profiles 
from  Figure  5-4,  first  locate  the  required  value  of 
fzQ/(0.185  u^) . The  geostrophic  wind  is  the  vector 
from  this  point  t.o  the  central  point  of  the  wind  spiral. 

Any  other  vector  drawn  from  this  point  to  any  point  on  the 
curve  represents  the  wind  velocity  at  the  height  indicated 
by  the  end  point. 
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FIGURE  5-4  hodograph  of  dimensionless  velocity  for  the  neutral 

BAROTROPIC  TURNING  LAYER  15-11 J. 


As  an  example,  suppose  zQ  = 0.01m,  f = 0.84  x 10"^  s 
and  u*Q  = 0.22  m/s,  the  geostrophic  wind  is  given  by 


where  W and  W are  determined  from  Equations  5.30  and 
yy 

5.31  with  h = 0.3u*/|f|.  Thus,  Wxg  = 5.11  m/s,  Wyg  = 1.63  m/s 

and  Wg  = 5.36  m/s.  The  wind  speed  and  direction  at  z = 162  m 

■>s  obtained  by  scaling  from  Figure  5-4  which  gives 

W = Scaled  Distance  (0.32-0.00002)  W /Scaled  Distance 
(center  - 0.00002)  9 

= (4.75)  (5.36  m/s) /5 . 16  = 4.93  m/s 

and  a = 7°  . 

Solutions  for  the  diabatie  barotrcpic  wind  profiles 
have  been  numerically  computed  [5-20,  21]  but  are  not 
expressible  in  a non-dimensional  form  such  as  Figure  5-4. 

The  experimental  results  of  Clarke  and  Hess  [5-16] 
provide  a method  of  estimating  the  wind  profile  for 
conditions  departing  from  those  of  neutral  stability. 

Figures  5-5  and  5-6  give  the  departure  of  the  nondimensional 

A 

wind  components  parallel,  (z) /u,Vo , and  perpendicular, 

A 

W (z)/u,v  , to  the  surface  wind  respectively  from  the 
reference  wind  defined  at  z = zf/u>VQ  = 0.15  as  a function 

A 

of  nondimensional  height,  z.  Using  the  values  of  A^(p) 

and  B x (u)  -where  \>=  ku*/|f|L,  given  in  Figure  5-7  and  5-8 

which  are  based  on  z = 0.15  and  Equations  5.30  and  5.31 

with  h =!f|/u,v0  the  values  of  Wxref  = Wx(z  =-  0.15)  and 

VJ  . = W t'z  = 0.15)  can  be  determined.  For  a given 
yref  y ' 
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stability  condition  y the  values  of  AWX  - (Wx(z  " °-15) 

- Wx(i)  )/u*Q  and  4Wy  - <W  <2  - 0.15)  - Wy(2))/u*c  can  be 

scaled  from  Figures  5-5  and  5-6  at  any  height  z up  to  a 
value  of  0.15. 

Figure  5-9  shows  a plot  of  the  wind  spiral  for  a 
stable  atmosphere  of  L*  11m  determined  in  the  above  manner. 

One  observes  that  the  turning  of  the  wind  begins  at 
approximately  20m  indicating  that  the  turning  layer  can 
become  very  significant  under  stable  nighttime  inversion. 

Figure  5-10  compares  the  log-linear  wind  profile  Equation  5*14 
and  Equation  5.11  with  the  horizontal  component  of  the 
turning  layer.  The  log-linear  law  begins  to  depart 
appreciably  at  20m. 

5.3.3  Baroclinic  Turning  Layer 

The  baroclinic  turning  layer  has  been  solved 
numerically  by  Blackadar  [5-22]  for  the  assumption  of 
a linearly  varying  geostrophic  wind. 

*G  - %o  + U 5'41 

where  VL,  is  the  surface  value  of  the  geostrophic  wind 
Go 

and  A is  the  vectorial  rate  of  increase  with  height.  The 
results  of  calculations  for  four  different  directions  of 
A are  shown  in  Figure  5-11.  The  curves  are  calculated 
for j A j = 42  x 105f,  |W  | = 84  x 105f,  zQ«  0.01m,  and  the 
indicated  heights  are  expressed  in  meters.  The  baroclinic 
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TO  THE  SURFACE 


EXPERIMENTAL  VALUES 


I MENTAL  VALUES 


FIGURE  5-9  CALCULATED  WIND  SPIRAL  FOR  A STABLE  ATMOSPHERE^ L~llM . 


FIGURE  5-11  NUMERICALLY  COMPUTED  BAROCLINIC  WIND  HODOGRAPHS  FOR 
FOUR  DIFFERENT  DIRECTIONS  OF  THE  THERMAL  WIND,  15-12]. 


TC^spa^y^TV-r 


wind  profile  can  be  scaled  from  these  curves  only  for 
the  specific  conditions  stated. 

Application  of  Equations  5.30,  5.31,  5.36  and  5.37 
requires  a different  scaling  velocity  than  WG<  Yamanda 
[5-15]  uses  a vertically  averaged  geostrophic  wind 


5.42 


Scaling  velocities  used  by  other  investigators  are  also 
reported  in  [5-15]. 


5 . 4 Boundary  Layer  Over  Non-Homogeneous  Terrain 
5.4.1  Surface  Roughness  Transition 

An  internal  boundary  layer  occurs  within  the  atmos- 
pheric boundary  layer  when  the  moving  air  mass  passes  over 
a change  in  surface  roughness.  Typically  the  interna] 
boundary  layer  will  develop  as  shown  in  Figure  5-12. 

When  a change  from  surface  roughness  to  surface 
roughness  zq2  occurs  at  location  x = 0,  a transition 
zone  grows  upward  and  spreads  outward  from  the  discontinuity 
in  surface  roughness.  Above  region  XI  shown  in  Figure  5-12 
the  wind  profile  is  undisturbed  from  the  upstream  profile 
which,  depending  on  stability  conditions,  is  described  by 
Equations  5.20  through  5.24.  Below  region  II,  the  wind 
speed  profile  may  be  taken  as  logarithmic 
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5.43 


, u*o2(x) 

x k 
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Further  from  the  discontinuity,  where  z/L  within  region  I 
becomes  large,  the  wind  speed  profile  may  more  adequately 
be  described  for  the  existing  diabatic  conditions  by  the 
appropriate  equations  given  in  Section  5.2.  Within 
region  II  the  wind  blends  gradually  from  the  velocity 
profile  in  region  I to  that  in  region  III.  At  large 
distances  from  either  side  of  the  discontinuity  in 
roughness,  the  shear  stress  at  the  surface  adjusts  to 
values  consistent  with  those  expected  of  the  flow  had  it 
developed  over  uniform  terrain. 

Numerous  mathematical  model  and  numerical  solutions 
of  the  internal  boundary  layer  have  been  proposed. 

Surveys  of  these  are  given  in  f5-23,  24,  25] . The 
experimental  data  to  confirm  these  models  are  somewhat 
limited,  however.  Reported  data  from  measurements  in 
the  laboratory  and  in  the  atmosphere  are  given  in 
1 5-26 , 27,  28,  29].  These  data  tend  to  corroborate 
the  mathematics,  but  are  generally  not  measured  over 
sufficient  distances  to  verify  the  character  of  the 
boundary  layer  at  elevations  where  encounters  with 
airplanes  occur.  In  the  absence  of  such  experimental 
data,  mathematical  models  are  used  to  extend  the  data  to 
heights  of  interest  for  this  report. 
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It  is  generally  agreed  both  experimentally  and 
analytically  that  the  height  of  the  internal  boundary 
layer,  essentially  6m  in  Figure  5-12,  grows  as: 


§^1-  e^jS-)" 

zo2  1 zo2 


5.44 


For  neutral  conditions,  values  of  n range  from  0.70  to  0.80 
and  values  for  c^  are  given  as:  0.75  - 0.03  In  (Z02/Z01^ 
by  Elliott  (5-30]  and  unity  by  Munro  and  Oke  [5-27], 


Rao  (5-31]  numerically  computes  for  unstable  flow  that  6 

Q GO  -I  q q 

grows  as  x for  L =-20m  and  as  x for  L =-2m.  For 
neutral  conditions  he  reports  n = 0.77. 


Figure  5-13  illustrates  the  approximate  growth  of 
the  internal  boundary  layer,  based  on  Equation  5.44 
relative  to  a 3°  glide  slope  for  some  typical  roughness 
changes  near  airports.  Wind  shear  may  be  expected  as 
aircraft  pass  through  the  internal  boundary  layer  and 
knowledge  of  the  shape  of  the  wind  profile  is  required  to 
assess  the  magnitude  of  this  shear. 

Solutions  of  the  two-dimensional  flow  field  over  a 
change  in  surface  roughness  by  numerical  solution  of  either 
the  boundary  layer  or  Navier-Stokes  equations  using 
various  turbulence  models  are  given  in  [5-12,  32,  and  34]. 
Typical  horizontal  velocity  profiles  at  various  stations 
in  the  windward  direction  from  the  surface  roughness  change 
are  shown  in  Figure  5-14,  and  a characteristic  vertical 
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Figure  5-13  COMPUTED  GROWTH  OF  THE  INTERNAL  BOUNDARY  LAYER 
RELATIVE  TO  A 3°  GLIDE  SLOPE  (DASHED  LINES)  FOR 
SOME  TYPICAL  ROUGHNESS  CHANGES  NEAR  AIRPORTS. 
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FIGURE  5-15  CHARACTERISTIC  VERTICAL  VELOCITY  FIELD  OVER  A 
CHANGE  IN  SURFACE  ROUGHNESS  [ 5“24l  . 


ev"a' 


velocity  field  is  shown  in  Figure  5-15.  The  velocity 
information  given  in  the  two  figures  are  for  a rough  to 
smooth  transition  of  In  z^/z c2=^*  Similar  results  are 
obtained  for  smooth  to  rough  transitions  and  some  results 
for  the  unstable  case  have  also  been  reported.  One  observes 
a downward  motion  of  the  air  as  it  passes  over  the  change 
of  surface  roughness.  This  suggests  an  unexpected  downdraft 
may  be  encountered  by  aircraft  approaching,  say  over  a 
forested  area  breaking  into  a cleared  airfield. 

The  application  of  the  above  numerical  solutions  to 

computer  simulation  of  airplane  flight  over  transitions 

in  surface  roughness  requires  either  on-line  solution  of 

the  complete  ser  of  equations,  which  is  extremely  costly 

in  computer  time,  or  a table  look-up  scheme  with  the  wind 

* 

field  data  stored  on  tapes  or  cards.  The  latter  approach 
is  readily  carried  out.  However,  for  quick  preliminary 
analyses,  a simple  mathematical  expression  for  the  wind 
profile  through  the  internal  boundary  layer  would  be 
valuable. 

Logan  and  Fichtl  (5-35 J have  recently  proposed  such 
a model  for  neutrally  stable  conditions  consisting  of  two 
bounding  logarithmic  layers  and  an  intermediate  velocity 
defect  layer.  Velocities  and  stress  distribution  functions 
within  region  II,  Figure  5-12,  which  meet  all  boundary 
and  matching  conditions  of  regions  I and  III  are  given  by: 
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u(z)  * "*(z,x)  F(n)  + ut 

vl  (z,x)  - (u*o2  - U lol)  G(n)  + U lol 


5.45 

5.46 


where  n»  z/6Q  and  6q  is  the  outer  boundary  of  region  II. 
The  velocity  u^  is  given  by: 


ui 

^ol 


5.47 


where  <z  <<$o  defines  a region  in  which  the  flow  has 
decelerated  relative  to  the  original  equilibrium 
logarithmic  profile: 


Wx  _ 1 •.  z 

“*01  " k ?ol 


5.48 


and  6.  <z<  defines  a region  in  which  the  flow  has 
x m ° 

accelerated  relative  to  the  logarithmic  velocity  profile: 


^ = I In  - 

u*o2  k zo2 


5.49 


By  appropriate  matching  techniques  the  relationships  given 
below  can  be  found: 


V {m  (1  + 5'50 

* - ’'yi  (r(ln  r-(l+M))  + em  ) 5.51 

z k x 


M = ln(zol/zo2) 


sm/zo2 


5.52 
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1 


5.53 


U*o2 


u 


*ol 


M 


ln(6m/2ol> 


+ 1 


u, 


u*ol  , 

T-  ln  > 


m 

ol 


5.54 


These  relationships  are  based  on  the  functions 

4(6  +6t-2z) 

GCn)  - erf  [ — ] 5.55 

where 

1 

erf<o  =/z? 

and 

F(n)  - p In 

The  expression  relating  6^,  6m  and  6q  is: 

Sm  = (<So  +{i>  12  5-50 

The  model  is  thus  complete  except  that  6^  must  be 
determined  from  experiment.  In  Reference  [5-35J  the  data 
of  Bradley  [5-26]  was  used  to  determine  6^.  The  authors 
of  [5-35]  point  out  that  6m  must  be  determined  separately 


e^2dS 


z 

6. 


m 


5.56 


5.57 
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for  each  analysis  from  experimental  data  having  tha  required 
value  of  20]/2q 2* 

To  circumvent  the  requirement  of  individual  experimental 
results  for  each  analysis,  it  is  proposed  that  simply 
be  related  to  x by 


V ‘°'8  2o2^ 


5.59 


(see  Equation  5.44) 

A linear  expansion  around  is  assumed  which  gives 


4o*  <x  + c2>s„, 


5.60 


c2Sm 


5.61 


The  solution  for  the  velocity  profile  within  the  internal 
boundary  layer  is  then  straight  forward.  With  C2X  0.5  the 
computed  velocity  profiles  using  Equations  5.52,  5.54, 

5.55  and  5.57  were  evaluated  and  are  compared  with  Bradley's 
data  in  Table  5-2. 


Table  5-2 


Comparison  of  Modified  Logan  and  Fichtll  5-35]  solution  with 
Bradleys  Data  [ 5-26] . 


x « 1220  cm 

MODEL 

BRADLEY 

z 

(cm) 

u(z)/uref 

u(z)/ur 

100 

0.96 

0.97 

60 

0.90 

0.895 

42 

0.87 

0.86 

x * 210  cm 

MODEL 

BRADLEY 

z 

(cm) 

uOO/uref 

u(z) /uref 

31.2 

0.79 

0.79 

20 

0.72 

0.72 

10 

0.67 

0.64 

The  agreement  is  seen  to  be  very  good.  Thus,  the  modified 
approach  to  the  solution  technique  of  Logan  and  Fichtl  [5-35] 
provides  an  easily  programmable  mathematical  model  for 
computing  internal  boundary  layer  wind  profiles.  Figure  5-16 
illustrates  computed  wind  profiles  at  various  stations 
downstream  of  a discontinuity  in  surface  roughness.  It 
appears  from  inspection  that  the  wind  shear  is  probably 
negligible,  however,  this  will  be  confirmed  by  introducing 
the  model  into  the  flight  dynamic  computer  code  in  a 
later  section. 
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COMPUTED  WIND  PROFILES  AT  VARIOUS  STATIONS 
DOWNSTREAM  OF  A DISCONTINUITY  IN  SURFACE  ROUGHNESS. 


5.4.2  Other  Surface  Variations 

The  preceeding  discussion  considers  only  winds  over 
changes  in  surface  roughness  where  the  wind  approaches 
over  a large  uniform  plane  and  re-establishes  itself  over 
an  equally  large  uniform  plane.  More  complex  boundary 
layers  may  develop  over  irregularly  spaced  patches  of 
roughness  such  as  illustrated  in  Figure  5-17,  see  [5-36]. 
These  types  of  wind  fields  may  best  be  analyzed  as  wind 
barriers  for  which  an  extensive  review  is  given  in  [5-37], 
When  the  surfaces  of  different  roughness  are  also  at 
different  temperature,  recirculation  can  occur  as  in  sea 
breezes  [5-38,  39,  40].  Sea  breezes  have  been'  reported 
to  cause  difficult  landing  conditions  [1-9].  Since 
many  airfields  are  located  near  the  sea  or  near  other 
large  bodies  of  water,  a model  of  sea  breezes  for  use  in 
flight  simulation  studies  should  be  developed. 

The  atmospheric  boundary  layer  over  "heat  islands", 
that  is,  surfaces  of  different  heat  flux,  temperature  and 
roughness  such  as  a large  metropolitan  area  surrounded  by 
flat  grass  land,  are  currently  being  studied  [5-24,  40,  41], 
In  view  of  the  fact  that  many  approach  and  departure  flight 
paths  are  over  expansive  cities,  wind  shear  over  heat 
islands  is  another  area  which  requires  investigation. 


6 . 0  SUMMARY 


6 . 1 Introduction 

Physical  models  of  wind  shear  environments  potentially 
encounterable  by  aircraft  in  the  terminal  area  have  been 
reviewed.  Numerous  mathematical  models  and  data  sources 
describing  wind  shear  conditions  as  reported  in  the  litera- 
ture are  surveyed.  Table  6-1  summarizes  the  result  of  the 
survey.  Sections  6.2,  6.3  and  6.4  discuss  the  tentative 
conclusions  and  recommendation  for  the  three  major  areas  of 
wind  shear;  fronts,  thunderstorms  and  neutral  and  stable 
boundary  layers,  respectively.  Section  6.5  describes  the 
turbulence  models  to  be  used  in  establishing  the  primitive 
turbulent  wind  shear  profiles  and  Section  6.6  discusses 
application  of  the  models  to  simulator  operations. 

6. 2 Fronts 

The  search  for  data  sources  from  which  mathematical 
models  of  wind  shear  can  be  established  for  synoptic  scale 
warm  and  cold  fronts  is  still  being  carried  out.  Data  from 
the  500  m tower  at  NSSL/NOAA  Laboratories  in  Norman, 

Oklahoma  have  been  measured  and  these  measurements  are 
currently  being  reduced  to  a format  similar  to  that  given  by 
Goff  (4-2]  for  thunderstorm  gust  fronts,  see  Figure  4-10. 

It  is  expected  that  these  data  will  be  available  by  7*pril, 
1977  and  they  will  be  tabulated  and  used  with  a computer 
lookup  routine  as  described  in  Section  6.3  for  thunderstorm 
gust  front  data. 
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SUMMARY  OF  LITERATURE  SURVEY 
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6.3  Thunderstorm  Gust  Front 


The  most  extensive  thunderstorm  gust  front  data  for 

formulating  detailed  wind  speed  profiles  are  those  measured 

by  Goff  [4-2]  in  the  form  of  cross  sections  of  lines  of 

constant  vertical,  W , longitudinal,  W , and  lateral,  W , 

z x y 

wind  speed  components  through  the  thunderstorm  parallel  to 
its  mean  motion.  Tabulating  these  data  on  a grid  and 
utilizing  a computer  lookup  routine  allows  the  wind  speed  to 
be  predicted  at  any  horizontal  and  vertical  position  called 
for  by  the  simulator  control  program. 

The  range  of  data  is  500  m in  vertical  and  3,321  m- 
11,685  m in  horizontal  extent,  depending  on  the  data  set* 
considered.  Twenty  data  sets  are  available  and  have  been 
punched  on  computer  cards  for  the  table-lookup  routine. 

Figure  6-1  illustrates  the  grid  system  superimposed  on  the 
wind  field.  Tabulated  values  of  wind  speed  for  the  20 
thunderstorm  gust  front  cases  are  given  in  the  appendix  in 
Tables  A-l  to  A-20. 

The  vertical  grid  spacing  is  the  same  for  all  tables, 
however,  the  horizontal  grid  spacing  varies  according  to  the 
length  of  record  available.  The  vertical  grid  spacing  DZ 
and  horizontal  grid  spacing  DX  are  listed  in  meters  at  the 
upper  right=hand  corner  of  the  table.  Each  data  set  has  11 
nodes  in  the  vertical  and  41  nodes  in  the  horizontal;  not  all 
nodes  are  shown  on  Figure  6-1.  Thus  for  example,  the  horizontal 
spatial  extent  of  the  record  in  Table  A-l  is  40  x 100  m = 4000  m 
whereas  in  Table  A-7  it  is  40  x 196  m = 7840  m. 
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FIGURE  6-1  VERTICAL  VELOCITY  CONTOUR  GIVEN  BY  f»OFF 

COMPARED  WITH  TABULATED  VALUES  FOR  COMPUTER 
LOOK-UP  GRID  SYSTEM, 


Figure  6-2  shows  the  flight  path  of  an  airplane 
characteristic  of  a DC-8  landing  with  fixed  controls  through 
the  thunderstorm.  Figure  6-3  gives  the  wind  speed  encountered 
by  the  airplane  during  descent.  The  wind  speeds  encountered 
are  seen  to  have  a very  complex  profile  indicating  that 
mathematical  formulations  of  gust  fronts  which  express  the 
wind  speed  only  as  a function  of  elevation,  such  as  Equations 
4.9,  4.10,  and  4.4  through  4.7  cannot  provide  realistic 
simulations. 

In  principle,  three-dimensional  wind  fields  are 
required  to  give  meaningful  simulations,  but  these  data  are 
not  available.  Thus,  simulation  of  flight  across  the 
direction  of  travel  of  the  storm  is  not  at  this  time 
possible.  Note  that  all  ci*oss  section  given  by  Goff  [4-2] 
are  in  planes  parallel  to  the  direction  of  travel  of  the 
storm  and  any  simulation  employing  these  data  assumes  the 
aircraft  is  flying  either  toward  or  away  from  the  storm 
along  the  direction  of  its  mean  motion.  Work  is  continuing 
to  provide  estimates  of  three-dimensional  gust  fronts  from 
Goff 1 s data  [4-2] . 

Developing  the  statistics  to  provide  risk  of  exceedance 
criteria  promises  considerable  difficulty.  The  approach 
envisioned  at  this  time  will  first  attempt  tc  establish  a 
common  point  on  all  20  storm  cases  such  as  the  location  of 
the  gust  front.  Second,  the  mean  wind  speed  for  all  20 
cases  will  be  computed  at  corresponding  grid  points  relative 
to  the  common  point.  Third,  the  standard  devi'i_ion  about 


100 


rlbUKL  b-Z  COMPUTERIZED  FLIGHT  PATH  OF  AIRCRAF' 
OF  D C 8 THROUGH  CASE  9 THUNDERSTO 


the  mean  with  the  assumption  of  normality  will  be  used  to 
estimate  the  probability  of  exceeding  a prescribed  wind 
speed  or  wind  shear  at  each  grid  point. 

6.4  Atmospheric  Boundary  Layer 

The  data  of  Clarke  and  Hess  [5-16]  from  Figures  5-5  and 

5-6  have  been  tabulated  in  Table  6-2  as  a function  of 

dimensionless  height  2 « fz/u * and  stability  parameter 

y m ku*/fL.  Coupled  with  a computer  lookup  scheme  the 

longitudinal,  W , and  lateral,  W , components  of  wind  speed 
x y 

can  be  determined  for  dimensionless  heights  of  0 < 2 < 0.15 
and  for  given  atmospheric  stability  conditions  of  -300  < 
y < 200.  The  range  in  physical  height,  z,  of  the  data  is 
dependent  on  the  latitude  through  the  Coriolis  parameter,  f, 
and  the  friction  velocity,  u*.  For  typical  values  of  u,(  - 
10~1,8~1,  z ranges  from  0 to  750  m. 

Wind  speed  profiles  generated  by  the  computer  lookup 
scheme  are  shown  in  Figures  6-4  and  6-5.  Note  that  the 
strongest  speed  shear  occurs  under  very  stable  conditions, 
y = 200,  but  the  strongest  directional  shear  occurs  under 
intermediately  stable  conditions,  y = 50. 

Actually,  the  strongest  directional  shear  occurs  under 
baroclinic  conditions.  A computer  code  for  computing  wind 
speed  profiles  under  baroclinic  conditions  where  the 
synoptic  pressure  gradients  vary  linearly  with  height  [5-22] 
is  programmed.  Realistic  magnitudes  for  the  rates  of  varia- 
tion, however,  are  being  justified  before  meaningful  wind 
speed  profiles  can  be  presented  for  simulation  work. 


103 


LONGITUDINAL  AND  LATERAL  WIND  COMPONENTS  AS  A FUNCTION 
OF  DIMENSIONLESS  HEIGHT  Z AND  STABILITY  u 
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Horizontal  wind  ahear  for  the  atmospheric  boundary 
layer  models  is  not  being  considered  in  this  program. 

Figure  5-16  suggests  that  horizontal  shear  resulting  from 
terrain  features  diffuse  rather  rapidly  and  should  not  cause 
significant  effect  above  elevations  of  approximately  25  m. 

Under  high  wind  conditions , however,  these  low  level  wind 
shears  due  to  vortices  or  recirculation  zones  near  buildings 
and  other  obstructions  to  the  surface  wind  may  create 
hazardous  flight  conditions.  This  is  an  area  where  a follow 
on  effort  to  establish  mathematical  models  of  hazardous 
wind  conditions  for  simulation  of  flight  near  building  is 
needed. 

Risk  of  exceedance  predictions  for  the  -atmospneric 
boundary  layer  require  establishing  the  probability  of  a 
given  value  of  y.  Figure  6-6  from  [5-16]  provides  some 
insight  in  this  regard.  The  figure  shows  a 40  day  average 
of  the  daily  variation  of  y with  two  standard  deviation 
error  bands.  Assuming  a Gaussian  distribution,  the 
probability  of  exceeding  a given  value  of  y can  be  estimated. 

This  curve  is,  of  course,  highly  sensitive  to  location  and 
further  research  is  required  to  establish  a more  general 
daily  variation  in  atmospheric  stability. 

6 . 5 Turbulence  Models 

The  primitive  turbulent  wind  speed  profiles  will  be 
generated  by  superimposing  a randomly  generated  turbulence 
signal  on  the  steady-state  wind  speed.  The  review  of  the 
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literature  pertaining  to  turbulence  model  is  not  included  in 
this  interim  report,  however,  the  reader  is  referred  to  [6-1 
through  7]  for  information  on  current  models. 

The  development  of  a realistic  turbulence  model  depends 
primarily  on  knowledge  of  the  appropriate  turbulence  energy 
spectra.  For  atmospheric  boundary  layers,  well  established 
turbulence  spectra  are  given  in  [5-2,  5,  6 and  6-8]  and 
turbulence  simulation  in  these  cases  is  straightforward. 

For  turbulence  simulation  in  thunderstorms,  spectra  are 
still  in  the  development  stage.  In  fact,  thunderstorm 
turbulence  is  neither  statistically  stationary  nor  homo- 
geneous and  superimposing  a turbulent  component  of  wind  on 
the  quasi-steady  state  values  given  in  Tables  A“1  through 
A -20  may  not  be  realistic.  However,  in  lieu  of  a more 
reliable  method,  this  approach  will  be  employed  until  the 
state  of  the  art  in  thunderstorm  turbulence  modelling  is 
improved. 

6 . 6 Flight  Simulator  Applications 

Mathematical  models  for  programming  simulators  to 
respond  to  fronts  are  not  ready  at  this  time.  Models  for 
steady-state  thunderstorm  gust  fronts,  as  described  in 
Section  6.3,  and  for  atmospheric  boundary  layers,  as 
described  in  Section  6.4,  are  available  as  computer  lookup 
routines  which  are  ready  for  immediate  application.  Addi- 
tional work  is  required  to  increase  the  sophistication  of 
the  models  to  include  turbulence  and  a statement  relative  to 
the  risk  of  encountering  a specific  wind  shear. 
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The  thunderstorm  gust  front  model  consists  of  20 
separate  data  sets.  The  suggested  application  of  these  data 
to  programs  for  evaluating  visual  displays,  for  testing 
avionics,  for  verifying  operational  procedures  or  for 
designing  control  systems  where  a pilot  poll  is  taken  as  the 
evaluation  criteria,  is  as  follows.  Each  pilot  should  "fly" 
as  many  of  the  wind  shear  cases  as  desired,  say,  10  to  20 
cases,  in  a prescribed  order.  The  pilots  would  be  unable  to 
learn  the  wind  shear  encountered  as  each  one  would  be 
different,  however,  the  experimental  sampling  would  be 
consistent  in  that  every  member  of  the  team  would  experience 
the  exact  same  sequence  of  wind  shears. 

A similar  procedure  can  be  carried  out  'with  the  atmo- 
spheric boundary  layer  model.  In  this  case,  a matrix  of  u 
values  would  be  established  prior  to  the  experiment  and  the 
simulator  would  be  programmed  to  accept  these  in  a specified 
sequence.  The  selection  of  values  for  \i  can  be  guided  by 
the  data  in  Figure  6-6. 

It  should  be  cautioned  that  the  present  data  decks  are 
not  likely  to  include  the  extreme  wind  shear  that  might  be 
encountered  once  in  a lifetime.  The  ability  to  simulate 
extremes  is  contingent  on  developing  the  risk  of  exceedance 
statistics  which  are  currently  under  investigation. 

Real  time  simulations  are  easily  achieved  with  the 
computer  lookup  routines  since  interpolation  between  tabu- 
lated values  requires  only  milliseconds  of  machine  time. 
Moreover,  data  card  decks  are  easily  incorporated  into 
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simulator  program  logic  [6-9,  10 J and  thus  facilitate 
application  to  existing  simulator  programs. 

6.7  Conclusion 

Efforts  are  continuing  to  develop  simulation  models  of 
detailed  wind  shear  profiles  in  large-scale  frontal  motions, 
to  establish  meaningful  turbulence  spectra  for  thunderstorms 
and  to  establish  risk  of  exceedance  criteria  based  on  valid 
statistics  of  extreme  wind  shears.  Models  of  quasi  steady- 
state  wind  shear  for  thunderstorm  gust  fronts  and  atmospheric 
boundary  layers  are  completed  and  ready  for  immediate 
application  in  wind  shear  hazard/flight  simulation  studies. 
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